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Fe 3 P 0403 forms a non-centrosymmetric lattice structure (space group R3m) comprising trian¬ 
gular motifs of Fe^^ coupled by strong antiferromagnetic interactions (|©cu"| > 900 K). Neutron 
diffraction from polycrystalline samples shows that strong frustration eventually gives way to an 
ordered helical incommensurate structure below Tn = 163 K, with the helical axis in the hexagonal 
ab plane and a modulation length of ~ 86 A. The magnetic structure consists of an unusual needle¬ 
like correlation volume that extends past 900 A along the hexagonal c-axis but is limited to ~ 70 
A in the ab plane, despite the three-dimensional nature of the magnetic sublattice topology. The 
small in-plane correlation length, which persists to at least T = TAr/40, indicates a robust blocking of 
long-range order of the helical magnetic structure, and therefore stable domain walls, or other defect 
spin textures, must be abundant in Fe 3 P 0403 . Temperature dependent neutron powder diffraction 
reveals small negative thermal expansion below Tn- No change in lattice symmetry is observed on 
cooling through Tn, as revealed by high resolution synchrotron X-ray diffraction. The previously 
reported reduced moment of the Fe^^ ions (aS'=5/2), ji ~ 4.2 /is, is confirmed here through magneti¬ 
zation studies of a magnetically diluted solid solution series of compounds, Fe( 3 _a,)Gaa,P 0403 , and 
is consistent with the refined magnetic moment from neutron diffraction 4.14(2) jib- We attribute 
the reduced moment to a modified spin density distribution arising from ligand charge transfer in 
this insulating oxide. 


I. INTRODUCTION 

Materials that display helical magnetic order, particu¬ 
larly those that are non-centrosymmetric, are of current 
interest for their potential technological uses and fun¬ 
damentally unique physical properties.® Materials with 
crystal structures lacking an inversion center, such as 
the B20-type metallic magnets MnSi,® Fe(i_a,)Coa,Si,® 
and FeGe,® as well as insulating Cu20Se03,® form he¬ 
lical (phase modulated) incommensurate magnetic or¬ 
dered phases which can be converted, by applying a small 
magnetic field, into a regular lattice of vortex-like topo¬ 
logical spin textures known as Skyrmions. Once created, 
Skyrmions are stable against dissipation due to their 
topological nature. Furthermore, they can be manipu¬ 
lated with electric currentJ^or thermal gradients,® mak¬ 
ing the transport of isolated Skyrmions promising as a 
new technology to realize high density and low dissipa¬ 
tion spintronics devices. Recent proposals for antiferro¬ 
magnetic (AFM) Skyrmions (those arising from an un¬ 
derlying AFM helical structure) have shown that these 
theoretically predicted objects may offer practical ad¬ 
vantages by avoiding the Skyrmion hall effect The 
role of defects in the “parent” helical structure appears 
to be essential in the generation of isolated Skyrmions, 
which can potentially be nucleated at sample edges,im 
at helical domain w all int er sect ions,and by combining 
helical domain walls 

In order to understand magnetic defects of helical 
magnets and their manipulation, it is therefore of interest 
to explore systems in which helical spin order is suscep¬ 
tible to defect formation. Helical order can be stabilized 
by the Dzyaloshinskii-Moryia (DM) interaction favor¬ 


ing non-coplanar spins, or alternatively through a bal¬ 
ance of competing interactions, i.e., frustration. In frus¬ 
trated systems, the ground state is usually selected from 
a plethora of closely competing phases; in principle this 
leads to flexibility in defect generation of helical struc¬ 
tures. For instance, substitution of non-magnetic atoms 
may allow for experimental control of defect formation, 
as interruptions in the connectivity of the magnetic sub¬ 
lattice can tip the delicate balance of frustrated interac¬ 
tions, potentially favoring commensurate structures or 
providing sites for defect nucleation. 

In this report, we present a reinvestigation of the 
triangle-based material Fe3P0403, known as the min¬ 
eral Grattarolaite.li^ We show that Fe3P0403 forms a 
high density of defects within an antiferromagnetic he¬ 
lical phase. Fe3P04 03 wa s first discovered and stud¬ 
ied via magnetization,^^^^^ specific heat, ^neu tron p ow- 
der diffraction,^^ and Mossbauer spectroscop}^ HUM sev¬ 
eral decades ago. Its non-centrosymmetric lattice (space 
group R 3 m) consists of triangular units of Fe^+ that are 
coplanar with the hexagonal ab plane and are linked to 
a PO4 group (Fig. a and b). Along the c axis, layers 
of triangular units arrange in a larger, triangular lat¬ 
tice pattern (Fig. [^c and d). However, the structure is 
not quasi-two dimensional: interplane (along the hexag¬ 
onal c axis) Fe-triangular units along the hexagonal c 
axis are separated by 3.13 A, while the intraplane dis¬ 
tance (in the ab plane) is 4.78 A. In effect, the three- 
dimensional magnetic sublattice is best described by a 
decorated “simple rhombohedral lattice,” i.e., a simple 
cube compressed along the body diagonal, which forms 
the “premedial lattice” composed of sites located at the 
center of the triangular subunits (Fig. [^e). 
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p re medial lattice 


FIG. 1. Schematic illustrations of the crystal structure of 
Fe 3 P 0403 , shown in the hexagonal setting of R3m. (a) The 
local coordination of Fe is a distorted trigonal bipyramid, 
forming FeOs clusters, (b) FeOs clusters are arranged in tri¬ 
angular subunits, with a PO 4 group linked to one face of the 
triangle, (c) The triangular units stack along c. (d) Sublattice 
structure of magnetic Fe^^ in Fe 3 P 0403 , viewed along the 
c-axis. Colored (grayscale) circles indicate the fractional co¬ 
ordinate (z) along the hexagonal c axis. Note that the overall 
phase along c is unconstrained by R3m. Here we choose pur¬ 
ple (dark gray, no line): z = 0.075, green (medium gray, left 
diagonal line): z = 0.409, orange (light gray, right diagonal 
line): z = 0.743. The two types of near-neighbor magnetic 
exchange interactions are shown in red (dark gray, Ji) and 
blue (medium gray, J 2 ) (e) Premedial lattice points (black 
spheres, located at the center of each triangular subunit) form 
a simple-rhombohedral lattice. 


Previous studies of Fe3P0403 reported a Neel transi¬ 
tion at T = 164.5 K via specific heaiP^land T = 173 T 5 K 
via magnetic susceptibility.^ Above this transition, the 
paramagnetic susceptibility was fit using two forms of a 
Curie-Weiss law, each of which gave unusual results for 
the effective moment, as will be discussed below. How¬ 
ever, it is clear that the interactions are strong (|0cw| ^ 
1000 K) and antiferromagnetic. The frustration param¬ 
eter, \&cw\/Tn^ is therefore greater than 6 and indi¬ 
cates significant frustration of the antiferromagnetic in¬ 
teractions. Below the transition, comparison of neutron 
powder diffraction at 4.2 K and at 200 K revealed the 
CO- development of broad fiat-topped peaks with a single 
sharp pea^ all of which were assumed to be magnetic 
in origin.l^ An AFM collinear commensurate magnetic 
structure was proposed, and was found to account for 
the central position of each diffraction peak and its in¬ 
tegrated intensity. The broadening of select peaks was 
attributed to an anisotropic correlation volume, but this 
was not studied further. 


In this contribution, we report new thermodynamic, 
synchrotron X-ray diffraction (SXRD), and neutron 
powder diffraction (NPD) measurements on Fe3P0403. 
The high momentum (Q) resolution of our NPD pat¬ 
terns allows us to rule out the previously proposed com¬ 
mensurate magnetic structure. We show that a helical 
incommensurate structure with a strongly anisotropic 
correlation volume best describes this magnetic neutron 
diffraction pattern. The model consists of magnetic do¬ 
mains restricted to 70 A in the hexagonal ab plane but 
unrestricted along the c axis; i.e., needle-like domains. 
Based on the magnetic structure, we estimate the ratio 
of the two nearest neighbor exchange couplings, allow¬ 
ing for analysis of the energy cost for simple domain 
walls. Although the microscopic mechanism stabilizing 
the high density of domain walls remains unclear, we 
demonstrate that the domain walls are expected to form 
perpendicular to the ab plane. 

The magnetic properties of the solid solution, 
Fe(3_a,)Gaa^P0403, in which non-magnetic Ga^+ re¬ 
places Fe^+, permits disambiguation of the size of the 
iron moment. Curie-Weiss analysis of high temperature 
susceptibility of the Fe(3_a,)Gaa;P0403 series provides a 
measure of the Fe^+ moment, which is consistent with 
magnetic model obtained from NPD. The moment is re¬ 
duced by 16 % from the expected value for S = 5 / 2 , 
pointing to covalency (ligand charge transfer) effects. Fi¬ 
nally, the temperature dependence of the chemical struc¬ 
ture is presented and shows some evidence for magneto- 
structural coupling, but without a symmetry-breaking 
structural phase transition occurring between 300 K and 
4.5 K. 


II. EXPERIMENTAL PROCEDURE 

Polycrystalline samples in the Fe(3_a,)Gaa,P0403 se¬ 
ries were synthesized by standard solid state meth¬ 
ods, following previous Fe203—FeP04 phase diagram 
reports.ll^ For the x = 0 compound, appropriate sto¬ 
ichiometric amounts of dried Fe203 and FeP04 were 
thoroughly mixed. Powders were pressed into pellets and 
sintered in alumina crucibles at 950 °C for 24 hours. Pel¬ 
lets were reground and heated multiple times at 1050 °C 
for 48 hours, until phase purity was maximized. Small 
amounts of Fe203 and FeP04 ( 1—4 % by mol) persisted 
even after three to four regrindings and reheatings. If 
residual FeP04 was detected (but only in the absence of 
Fe2 03), it was removed by overnight suspension in dilute 
HCl (2 M), followed by rinsing with deionized water and 
drying in air. 

Phase identity and cell parameters were determined 
by refining powder X-ray diffraction patterns from a 
Scintag Advanced Diffraction System (Cu-Ka radiation, 
A = 1.541 A) using the Rietveld method as implemented 
by GSAS/EXPGUL^ The remaining members of the 
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Fe( 3 _^)Ga^P 0403 series {x = 0.5,1.0,1.5,2.0,2.5,2.7) 
were prepared and analyzed similarly with stoichiomet¬ 
ric amounts of dry Fe 203 , Ga 203 , FeP 04 , and GaP 04 . 
Samples of high Ga—content {x = 2.0, 2.5, 2.7) required 
more heatings to achieve maximal phase purity, and the 
Fe( 3 _a,)Gaa,P 0403 type structure was not the dominant 
phase at x = 2.9 indicating a solubility limit between 
X = 2.8 and 2.9 (See Appendix [A|). None of the solid- 
solution compounds were treated with HCl. 


The unit cell parameters, a and c, for 
Fe( 3 _a,)Gaa,P 0403 were extracted from Rietveld 
refinements of powder X-ray diffraction patterns. These 
unit cell lengths decrease with increasing Ga(III) (ionic 
radius in 5-fold coordination= 55.0 pmp^ substitu¬ 
tion for Fe(III) (ionic radius in 5-fold coordination 
= 58.0 pm)^ (Fig. 11). As expected from Vegard’s law, 
the decrease in lattice parameters is linear as a function 
of X. 


Magnetization and specific heat measurements were 
performed using a Quantum Design Inc. Physical 
Properties Measurement System. The magnetization of 
Fe 3 P 0403 and the gallium-diluted Fe( 3 _a,)Gaa,P 0403 
series were measured from 1.8 K to 700 K under zero- 
field cooled (ZFC) and field-cooled (FC) conditions at 
/ioH = 1 T. Specific heat measurements were carried out 
using the semi-adiabatic heat-pulse method at tempera¬ 
tures from 1.8 K to 300 K. 



7-(K) 


FIG. 2. Cp/T vs. T (left axis) and M/H vs. T (H = 10000 Oe, 
right axis) for poly crystalline Fe 3 P 0403 . The data show a 
Neel transition at Tn = 163 K. The field cooled (FC) vs. zero 
field cooled (ZFC) magnetization measurements are nearly 
identical over the full temperature range. 


III. RESULTS 

A. Thermodynamic measurements 

1 . Specific Heat 


Neutron powder diffraction (NPD) experiments were 
carried out at T = [4.5, 60,100,120,160, 220, 295] K on 
a polycrystalline sample of Fe 3 P 0403 (m =1.8 g) using 
the BTl powder diffractometer at the NIST Center for 
Neutron Research. A constant wavelength of 2.0799 A 
with 15’ collimation was useT providing a resolution of 
SQIQ = 0.013 at Q = 1.0 polycrystalline sample 

of Fe 3 P 0403 (m =3.3 g) was studied at T = 295 K and 
100 K using the POWGEN time-of-fiight neutron diffrac¬ 
tometer at the Spallation Neutron Source, Oak Ridge 
National Laboratory. The long-wavelength frame (cen¬ 
tral wavelength = 2.665 A) was used for magnetic struc¬ 
ture refinements, providing a resolution of 5QIQ = 0.005 
at Q = 1.0 A. Synchrotron X-ray Diffraction (SXRD) 
data were collected on Fe 3 P 0403 using the diffractome¬ 
ter on the 11-BM-B beam line at the Advanced Photon 
Source, Argonne National Laboratory at temperatures 
of T = 100 K and 295 K using an incident wavelength 
of 0.414 A.l^ To account for X-ray absorption of the 
sample, a correction of /iR = 0.40, calculated from the 
sample composition and radius, was applied to the full 
pattern analysis 


The specific heat of Fe 3 P 0403 was originally reported 
in Ref. [T71 An anomaly was observed near 163 K, cor¬ 
responding to the magnetic transition temperature. We 
have confirmed this specific heat, shown in Fig.[^ In con¬ 
trast to Shi et al, we do not attempt to isolate the mag¬ 
netic contribution to the specific heat. Due to the pres¬ 
ence of strong magnetic correlations above Tat, evidenced 
by both inverse susceptibility (Fig.[^ and NPD (Fig.[^, 
estimating the lattice contribution is not accurate and 
appears to lead to spurious results such as the unphysi- 
cally small change in entropy reported in Ref. [iTl which 
is likely amplified over the broad temperature range of 
the anomaly in the specific heat. Our attempts to syn¬ 
thesize the X = 3 end member of Fe( 3 _a^)Gaa,P 0403 , 
Ga 3 P 0403 , for use as a lattice analog were not success¬ 
ful. 


2. Low Temperature Susceptibilty 

The anomaly observed in the specific heat of 
Fe 3 P 0403 coincides with an inflection point in its 
sharply decreasing magnetization, confirming the forma- 
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FIG. 3. (a) Susceptibility measurements (M/H) 

made on representative members of the solid solution 
Fe( 3 _a,)Gaa^P 0403 with an applied field of 10000 Oe. At 
X = 0, i.e., Fe 3 P 0403 , the FG and ZFG curves do not fully 
diverge, but show a small splitting around Tn- At higher 
X, the FG vs. ZFG curves split at a temperature which 
decreases with increasing x, as might be expected for a 
spin glass. The temperature of the splitting is referred to as 
Tp. (b) Tn and Tp as a function of gallium content (x) in 
Fe( 3 _a,)Gaa;P 0403 . Solid line is a guide to the eye. Dotted 
lines show possible trends of Tw and Tp as a function of x, 
but are not constrained by the data. The blue (gray) region 
indicates composition range for which no distinguishing 
features are observed above T = 1.8 K (based on data taken 
at the two compositions marked by open blue squares). 


tion of an ordered antiferromagnetic structure at 163 K 
(Fig. [^. The ZFC and FC susceptibility at /ioH = 1 T 
do not diverge appreciably, though there is a small gap 
between the two which develops above T/v, is maximal 
at T/v, and takes on a nearly constant value down to 1.8 
K. There is an upturn in the magnetization near 1.8 K 
which likely arises from “orphan spins” at lattice defects. 

Members of the solid solution Fe( 3 _a,)Gaa,P 0403 show 
a dramatic change in the the behavior of their low tem¬ 
perature susceptibility compared to Fe 3 P 0403 , as shown 
in Fig. for members x = 0.5 and x = 1.0. At x = 0.5, 
the AFM transition appears to be totally suppressed, 
but a transition which resembles freezing (due to the 
presence of a ZFC / FC splitting) occurs at Tp = 40 K. 
This feature persists at higher x but Tp decreases; at Ga 
concentrations higher than x = 1.5, any such features. 



FIG. 4. (a) and (b) Results of Gurie-Weiss analysis of the 
inverse magnetic susceptibility of Fe( 3 _a;)GaccP 0403 . Fitting 
was carried out in the temperature range of 300 K to 700 K. 
For X > 1.0 (red diamonds), the analysis produces unreliable 
results due to non-linearity in At higher x, the magnetic 
dilution the inverse susceptibility data become linear over 
300-700 K to allow reliable application of the Gurie-Weiss re¬ 
lation. As the range of compositions increases x > 1.0, the 
effective moment tends towards peff = 5.1 fiB, less than the 
value expected for S = 5/2 (5.9). (c) and (d) Inverse sus¬ 
ceptibility at H = 10000 Oe for x = 0 (c) and x = 2.5 (d). 
Orange lines are fits to the Gurie-Weiss law. 


if they exist, are outside of the temperature range of 
the measurement (i.e., Tp < 1.8 K). Fig.j^b) illustrates 
the magnetic phase evolution as a function of x in the 
Fe( 3 _a,)Ga^P 0403 system. 


3. High Temperature Suseeptihility 

The magnetic susceptibility in Fe 3 P 0403 and 
Fe( 3 _a,)Gaa,P 0403 measured from 300 K to 900 K were 
analyzed with the Curie-Weiss relation (the susceptibil¬ 
ity for Fe 3 P 0403 over the full measured temperature 
range, 1.8 K to 900 K, is shown in Appendix 0. Con¬ 
stants C and Qcw were extracted rrable[I|) by fits to the 
inverse Curie-Weiss relation (Eqn. [r via the linear least- 
squares method from 300—700 K. The effective paramag¬ 
netic moment, Peff = \/8 C, was subsequently calculated 
using the extracted C value. 


1 _ (F - Qcw) 

C 


( 1 ) 


We find that samples of Fe( 3 _a,)Gaa,P 0403 from x = 
0.0 to X = 1.0 display nonlinear inverse susceptibilities 
(Fig.0 and are unfit for Curie-Weiss analysis. Dilution 
with nonmagnetic Ga above 50% substitution from x = 
1.5 to X = 2.7 produces a paramagnetic response that is 
well described by the Curie-Weiss relation. The magnetic 
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TABLE L Results from magnetic susceptibility fits to 
the inverse Curie-Weiss equation 0 for Fe 3 P 0403 and 
Fe( 3 _a;)Gaa;P 0403 from 300 — TOOK. The effective moment 
is calculated as peff = VSC, where C is the inverse slope 
of the fit in CCS units of emU'K“^-molFe“^-Oe“^. Negative 
values of Gcw indicate an overall antiferromagnetic interac¬ 
tion. 


X 

C 

Peff 

(M-b) 

Ocw 

(K) 

1.5 

4.16 

5.77 

-545 

2.0 

3.72 

5.45 

-292 

2.5 

3.54 

5.32 

-97 


behavior of Fe atoms can be described by the extracted 
constants {Qcw^ C) and the relations C = (peff/8)^: 
where Peff = + 1) (Tabled. As x increases, the 

magnitude of &cw decreases, reflecting the decrease in 
the antiferromagnetic Fe-Fe interactions. The Fe atom 
effective moment trends to ^ 5.1 /i^ as &cw falls below 
the lowest temperature used for fitting (Fig. [^. Typical 
spin-only Fe^+ {S = 5/2) produces Peff = 5.9/iB- For the 
members of the solid solution that are reliably described 
by the Curie-Weiss relation {x > 1.5), Peff is noticeably 
lower than the spin-only effective moment. 


B. Synchrotron X-ray Diffraction on Fe 3 P 0403 

Structural refinements of SXRD data taken above (295 
K) and below (100 K) the magnetic transition temper¬ 
ature were performed using the published crystal struc¬ 
ture as a starting configuration.^ No change in sym¬ 
metry was observed on cooling through T/y. A 3.56wt% 
Fe 203 impurity was included in the Rietveld refinements. 
The SXRD data and refinements are shown in Fig. 

The crystal structure of Fe 3 P 0403 consists of three 
crystallographic sites for oxygen and unique sites for Fe 
and P. Note that in the hexagonal setting of R3m the 
overall fractional 2 ; coordinate is not fixed by any of the 
special positions occupied by atoms in Fe 3 P 0403 . We 
therefore chose an arbitrary fixed value of 2 : = 0.7425 for 
Fe. Table [ll] shows the refined positional and unit cell 
parameters, which display little variation from 295 K to 
100 K, and the isotropic displacement parameters (Uiso) 
decrease upon cooling to 100 K, as expected. 


C. Neutron Diffraction 

1. Magnetic Structure 

Fe 3 P 0403 develops magnetic reflections upon cooling 
through Tat = 163 K. Sharp (resolution limited) mag¬ 


netic peaks coexist with broad, flat-topped peaks having 
widths approximately equal to eight (twenty) times the 
resolution of the BTl (POWGEN) measurement. The 
magnetic nature of these features is confirmed by com¬ 
parison to the low temperature SXRD, which shows no 
extra peaks on cooling through the magnetic transition 
temperature (Fig. |^. 

The center of all magnetic features can be approxi¬ 
mately indexed to a = (0,0,1.5) ordering wavevec- 
tor, where the “h” subscript refers to the hexagonal 
setting of the unit cell (in the rhombohedral setting, 
kr = (0.5,0.5,0.5)). Such a wavevector produces a tt 
phase shift of the magnetic structure for each triangu¬ 
lar layer along the hexagonal c axis. In order to repro¬ 
duce the breadth of the flat-topped peaks, a modula¬ 
tion in the ab plane must be added, i.e., = (^a, ^ 6 , 0 ). 

This splits reflections at zone centers that have nonzero 
components in the ab plane, but does not split reflec¬ 
tions having a purely c-axis component. For instance, 
kh = (0, 0, 1.5) becomes (^^,^6,1.5) and therefore 
only slightly shifts in \Q\ for small 6a and 6i). In con¬ 
trast, reflections such as (1, 0, 0.5) and its 5 other 
symmetry-equivalent peaks become, e.g. (1 + (5^, (5^, 0.5), 
{Sa, 1 + 6b, 0.5), (1 + 6a, -1 + 6b, 0.5), etc., generating a 
series of six reflections that are closely-spaced in Q. 


Inspection of the high resolution POWGEN data re¬ 
veals that the shape of the broad peaks cannot be repro¬ 
duced using a sum of closely spaced resolution-limited 
peaks (Fig. [^. In fact, even a nearly continuous distri¬ 
bution of peaks, which could hypothetically result from a 
continuous degeneracy of in-plane modulation wavevec- 
tors, cannot account for the shape of these features (see 
the supplemental peakshape analysis in Fig. 12). This is 
evident from the slope of the edges of the broad features. 
In order to account for the peak shape, some broadening 
of in-plane peaks must be applied. This indicates finite 
sized domains in the a and b directions. In order to repro¬ 
duce the peak shape in the POWGEN data, we applied 
Gaussian broadening of FWHM = 0.03 to the low 
angle reflections generated by a helical structure with 
kh = (0.0259, —0.0902,1.5) (the same magnetic struc¬ 
ture refined from the BTl data, discussed below). This 
broadening is a factor of 6 greater than the instrumen¬ 
tal resolution, and corresponds to an in-plane correla¬ 
tion length of only ^ 70 A, while the helical modulation 
length is 86 A. Importantly, however, the c axis correla¬ 
tions remain long ranged, as indicated by the resolution- 
limited {6a, 6b, 1.5) peak, implying correlations longer 
than 900 A. Therefore, the diffraction pattern strongly 
suggests needle-like magnetic domains, with the longest 
correlations along the hexagonal c-axis, along which un¬ 
modulated antiferromagnetic correlations form. 


The shape of broad peaks does not change on cooling 
to at least 4.5 K (3% of T/v) (Fig.[^. This indicates that 
the short range magnetic domains are either equilibrium 
configurations or are kinetically inhibited from further 
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TABLE IL Crystal structure parameters of Fe 3 P 0403 (space group R3m) obtained from Rietveld analysis of 11—BM data 
(SXRD) at 100 K and 295 K. The z—position of the Fe-atom is held constant when rehning atom positions (see text). Uiso is 
shown in units of 10“^ A^. 


Unit Cell a 

Parameter c 


Atom 


Wyckoff 

Position 


X 


100 K 
7.999(1) A 
6.852(1) A 


y z 


Uu 


X 


295 K 
8.004(1) A 
6.860(1) A 

y z Uiso 


Fe 

9b 

0.7967(1) 

-0.7967(1) 

0.7425 

1.63(2) 

0.7970(1) 

-0.7970(1) 

0.7425 

4.83(3) 

P 

3a 

0 

0 

0.0015(1) 

1.79(12) 

0 

0 

0 .0002(1) 

5.21(13) 

01 

3a 

0 

0 

0.2309(2) 

2.27(30) 

0 

0 

0.2300(2) 

4.09(33) 

02 

9b 

0.5397(1) 

-0.5397(1) 

0.8503(1) 

3.32(18) 

0.5406(1) 

-0.5406(1) 

0.8482(1) 

6.79(21) 

03 

9b 

0.5620(1) 

-0.5620(1) 

0.2723(1) 

6 .01(21) 

0.5612(1) 

-0.5612(1) 

0.2712(2) 

14.37(27) 




FIG. 5. High resolution X-ray diffraction Rietveld pattern refinements of Fe 3 P 0403 at T = 100 K (left) and T = 295 K (right) 
collected at the Advanced Photon Source 11-BM-B beam line using an incident wavelength of 0.414 A. Observed data (x) 
are plotted with the Rietveld refinement profile (orange line). A 3.56 wt% Fe 203 impurity is present in the sample (bottom 
ticks). No change in lattice symmetry is observed between 100 K and 295 K. 


growth. Furthermore, the relevant wavevectors for the 
modulation do not change as a function of temperature. 
No higher harmonics are observed at any temperature, 
which is consistent with the absence of the “squaring up” 
of a longitudinally modulated structure.!^ The possibil¬ 
ity of a square-wave structure that develops immediately 
upon cooling through T/v was also investigated and found 
to not reproduce the peak shape (Fig.jl^. These consid¬ 
erations imply that the underlying magnetic structure is 
helical or conical with all ordered moments being of the 
same size. 

In order to further refine the magnetic structure of 
Fe 3 P 0403 , we used an ordering wavevector of the form 
((5a, (5^, 1.5) in the Fullprof Rietveld refinement pro¬ 


gram. Fullprof allows the refinement of anisotropic size 
broadening term with Lorentzian shape, accounting for 
either needle-like or platelet correlation volumes, for 
constant-wavelength NPD data.^^ We refined the Loren- 
tizan broadening term for a needle-like correlation vol¬ 
ume with long axis parallel to c, using the (constant- 
wavelength) BTl data at T = 4.5 K. 

The determination of the relative spin orientations 
proves difficult based solely on the powder diffraction 
data due to domain and powder averaging. However, 
certain possibilities can be ruled out. For instance, in or¬ 
der to generate the ^^,1.5) peak (the observed sharp 
magnetic peak), there must be a net moment in the ab 
plane in the structural unit cell. This rules out helical 
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FIG. 6. (a) Neutron powder diffraction from the POWGEN instrument at T = 10 K (< Tn) and at T = 300 K (> Tat). 
Broad magnetic peaks coexisting with a sharp peak appear below Tn- The broad features span ~ 18 times the instrumental 
resolution, while the sharp magnetic peak is resolution-limited. The centers of all magnetic features can be indexed to a = 
(0, 0, 3/2) ordering wavevector. (b) The width of the broad peaks can be partially accounted for by allowing an incommensurate 
modulation in the ab plane, kh = ((5a, 3/2). Gaussian broadening is applied to the lowest angle feature, which corresponds 

to several reffections that are nearly within the a*b* plane; this additional broadening is necessary to reproduce the slope of the 
edges of the peak, implying short-range order with Gaussian correlations in the ab plane with a correla tion len gth of ^ab = 70 
A. Note that the free ion Fe^^ magnetic form factor is not included in this comparison (see discussion in III G l| . Including the 
form factor does not affect the necessary broadening (see Fig. [^and related discussion). Inset: the {Sa,Sb) mociulation used to 
generate the (1 + (5a, (5?,, l/2)-type peaks shown in the main panel is represented by a white dot on a Toaap showing that the 
best fits are obtained with a modulation wavevector of modulus |(5| = 0.073 (blue circle). 


variations of 120°—type order on the triangular units. 
The best fits can be obtained by starting with a collinear 
magnetic structure, with moments aligned on each tri¬ 
angle and a tt phase shift from layer to layer (i.e., a com¬ 
mensurate (0,0,1.5) parent structure), and then allowing 
an a5-plane modulation to produce a helical structure. 
The best fit is obtained when the helical plane direction, 
n, lies in the ab plane; for example, the moments can 
be confined to the ac plane as in the structure refine¬ 
ment shown in Fig. Other types of phase-modulated 
structures cannot be ruled out, such as an antiferromag¬ 
netic conical structure with the cone axis along a and the 
opening angle being ^ 70°, though this type of structure 
would be less favored by Heisenberg exchange (Appendix 
We therefore choose the simplest structure, a helix 
as shown in Fig. for the purposes of further discus¬ 
sion. In all of these types of structures (with the cone or 
helical axis in the ab plane), the moment size refines to 
^ 4.1 to 4.3 liB^ which is significantly lower than the 5 
fiB expected from the S=b/2 state of Fe^+, but is con¬ 
sistent with a Curie-Weiss analysis of the magnet ically 
dilute series Fe( 3 _a,)Gaa,P 0403 (see Section III A3). 

In order to determine the direction of the modulation 
vector in the ab plane, S, we systematically investigated 
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a range of {Sa^Siy) pairs by refining the following: the 
two angular parameters, 0 and 0, specifying the direc¬ 
tion of the helical plane (n) relative to the c and a 
axes, the magnetic moment, and the anisotropic (needle¬ 
like) Lorentzian anisotropic broadening, using the Full- 
prof refinement software in the helical magnetic structure 
mode. The fitting range was restricted to Q = 0.27 to 
3.2 A“^ in the BTl data at T = 4.5 K. The resulting x^ 
values for the fits are shown in the plane in the 

inset of Fig. which shows that only the magnitude of 
the modulation wavevector can be constrained based on 
this fit. The blue line shows the modulus of best fit of 
\6\ = 0.087 ± 0.009 r.l.u.l^This shows that a particular 
modulation direction cannot be distinguished based on 
these data. 


Based on simple near neighbor interactions, we argue 
in Section |IV B| that there may actually be no preferred 
direction for the modulation vector in the system, i.e., 
a single crystal diffraction pattern would hypothetically 
show incommensurate rings of scattering. On the line of 
constant modulus which produces the lowest values, 
0 does not deviate from 90° by more than 4°, confirming 
that the helical plane contains the c axis, while the mag¬ 
netic moment does not deviate from 4.1 (ib by more than 
























0.2 and the Lorentzian broadening term varies from 
6.2 to 10.7, corresponding to correlation lengths of 45 to 
85 A. Visual inspection of various choices for 5 along this 
line of “best fit” favors a fit with 5 = (0.028,-0.097), 
shown in Fig.[^ and here the broadening term is 6.9, cor¬ 
responding to ^ah r\j 70 A, in agreement with our peak 
shape analysis of the POWGEN data using same mag¬ 
netic structure (Fig. |^. 

Finally, we note that the lowest angle broad peak is ob¬ 
served to be rather symmetric in Q, whereas the model, 
which includes the spherical form factor for Fe^+ as pa¬ 
rameterized in Ref. na does not fully reproduce this. We 
made many attempts to find a magnetic structure that 
would produce the symmetric shape of the lowest an¬ 
gle feature, while retaining other important features, but 
these attempts were not fruitful. We note that removing 
the form factor allows a better fit of the lowest angle 
feature (see Fig. [^. Although it is not physical to com¬ 
pletely neglect the form factor, a flatter form factor near 
this low angle peak can arise from covalent spin density. 
Similar effects have been studied in detail for other tran¬ 
sition metal ions; for instance, Ni^+ is known to have a 
more contracted spin density in the compound NiO com¬ 
pared to the free ion, leading to a more gradual decrease 
in the form factor.!^ Given the mounting evidence for 
liga nd cha rge (and spin) transfer in Fe3P0403(see Sec¬ 
tion IV A), the moment reduction and poor agreement 


with conventional magnetic form factors are therefore 
potentially related. Further measurements on single crys¬ 
tal samples will be required to understand details of the 
magnetic form factor, the domain anisotropy, and the 
spin orientations. 


2. Temperature Dependenee from Neutron Powder 
Diffraetion 

The thermal dependence of the basic structural pa¬ 
rameters were determined using the BTl data, by re¬ 
fining a nuclear and magnetic phase at each tempera¬ 
ture (taking the planar helix structure shown in Fig. 
as the magnetic structure). The lattice parameters (unit 
cell lengths a and c), cell volume, and ordered magnetic 
moment are shown as function of temperature in Fig. 
[^b). The a unit cell length and the cell volume vary 
non-monotonically below Tjy = 163 K; a begins to in¬ 
crease below ^ 100 K, exhibiting a minimum just below 
T/v- Meanwhile, c does not decrease enough to compen¬ 
sate the increase in a; a corresponding negative thermal 
expansion is observed below 100 K. These trends may 
be correlated with the onset of magnetic order, but the 
current measurements are not conclusive on this point. 
The other structural parameters, such as anisotropic dis¬ 
placement parameters, cannot be refined accurately from 
this data, as the experiment was optimized for magnetic 
(low Q) scattering. 


IV. DISCUSSION 


A. Curie-Weiss analysis and reduced Fe^^ moment 


Previous Curie-Weiss analysis of the magnetic suscep¬ 
tibility in Fe3P0403 found effective magnetic moments 
ofpeff = 6.45/i^BEland 4.3/iBl^and Curie-Weiss temper¬ 
atures of-1707 kP and -1087 kP With five electrons 
in five 3(i-orbitals, Fe^+ most often takes on a spin-only 
5'=5/2 angular momentum, with an expected effective 
moment of Peff = 5.9/iB- Therefore the magnetic mo¬ 
ments determined from the previous Curie-Weiss anal¬ 
ysis are anomalous and mutually inconsistent. Exami¬ 
nation of the inverse susceptibility in the temperature 
ranges used for fitting reveals why this is the case; the 
curve is non-linear up to at least 900 K. The obtained 
Curie-Weiss temperatures, Qcw >1000 K, give a hint 
as to why this may be: the temperature range used for 
fitting is exceeded by the deduced mean field interac¬ 
tion strength, invalidating the mean field approximation. 
Strong magnetic correlations above T/v are also revealed 
from the NPD patterns (Fig. [^; significant diffuse scat¬ 
tering around the (0, 0, 1.5) position is evident above 
160 K. This is expected in frustrated magnets. 

By this reasoning, the true mean field interaction 
strength, normally encoded by Qcwi is presently un¬ 
specified, but must be greater than ^ 900 K. The 
effective moment can still be deduced by diluting 
the magnetic lattice, as we have done using the 
Fe(3_a,)Gaa,P0403 series (Fig.0. The Curie-Weiss anal¬ 
ysis of that series shows that^crv decreases with in¬ 
creasing X (as expected), reaching a value that falls 
within our fitting range (300 K - 700 K) around x = 1.5. 
The effective moment per Fe^+ approaches Peff ^ 5.1 i^b 
at these higher x values. 

What could be the cause of this moment reduction? 
Orbital contributions are not expected in the half filled 
shell of Fe3+, so the ^f-factor is not likely to be sig¬ 
nificantly reduced. The most likely possibility is lig¬ 
and charge transfer, i.e., 2p(0^“) ^ 3d(Fe^+), which 
reduces the percentage of unpaired spins on Fe. Such 
effects are not uncommon and have been proposed 
for other Fe^+—containing compounds such as (a-Fe2 03 
(Hematite),!^ and Ca3Fe2Ge30i2!^ Indeed, such cova¬ 
lency effects were previously proposed for Fe3P0403.tI^ 

The reduced effective moment is consistent with a g- 
factor of 2 and a reduced total spin value of S' = 2.1 per 
Fe. This leads to a microscopic moment of /r = 4.2/iB, 
which is consistent with our refined value of the ordered 
magnetic moment from NPD (4.14(2) /i^), recalling that 
quantum fluctuations in antiferromagnets tend to further 
reduce the ordered moment. Thus, we conclude that the 
moment associated with Fe in Fe3P0403 is reduced by 
16% compared to the expected /i = 5/^^ for an isolated 
Fe^+ ion. 
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FIG. 7. Temperature dependent results of NPD from Fe 3 P 0403 taken on the BTl instrument, (a) Evolution of the diffraction 
profile near the lowest angle magnetic peaks, (b) Magnetic moment, vs. T from the magnetic structure rehnement shown 
in Fig. Solid line is a guide to the eye. (c) Unit cell volume in the hexagonal setting; a small negative thermal expansion 
is observed below 100 K. d) and e) Temperature dependence of the unit cell lengths a and c, respectively. The cell length, a, 
decreases on cooling but begins to increase below Tn- Errorbars are smaller than markers in b) through e). 
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EIG. 8. (a) Results of Rietveld refinement of the BTl neutron 
powder diffraction data at 4.5 K. The nuclear (top ticks) and 
magnetic (bottom ticks) phases are shown. Broadening of in¬ 
plane magnetic peaks was applied (see text). The difference 
curve is shown in the bottom box (blue line), (b) Illustration 
of of the refined helical magnetic structure, showing 10 unit 
cells superimposed along b. The highlighted arrows show the 
moment directions in the nearest unit cells. Thick black lines 
connecting atoms are Ji bonds (i.e., the edges of the trian¬ 
gular units), (c) View of the magnetic structure along the a 
axis, which is perpendicular to the helical plane. Note in b) 
and c) how moment directions are reversed in neighboring 
unit cells along the c-axis. 


B. Microscopic Origin of Helimagnetism 

Two common sources for helimagnetism are 1) the 
anti-symmetric Dzyaloshinskii-Moriya (DM) interaction, 
and 2) frustration arising from a competition between Ji 
and J 2 . We first discuss the effect of DM interactions in 
Fe3P0403. 

In Fe3P0403, DM interactions are allowed for each 
pair of Fe ions. Here we consider only nearest neigh¬ 
bors. A mirror plane bisects each nearest neighbor bond. 
AccordiM to symmetry arguments first articulated by 
Moriya,!^ this leads to a DM vector constrained to lie in 
the plane perpendicular to the Ji bond. The triangular 
geometry dictates that any component of the DM vector 
lying in the ab plane for a given near neighbor bond will 
then be cancelled by the remaining two, leaving only the 
possibility of a net DM vector along the c-axis. There¬ 
fore the nearest neighbor DM interaction could only fa¬ 
vor spin structures having a net pairwise cross product 
pointing along ±(0,0,1). We made many attempts to 
find a suitable spin structure based on this constraint, 
but neither a helical nor conical structure with a net 
cross product along c produces a satisfactory fit to the 
data. In fact, the best fits are obtained from configura¬ 
tions with a net cross product perpendicular to c. 

The most likely cause for the observed helimagnetism 
is therefore the competition between AFM Ji and J2. 
The type of magnetic order that develops in Fe3P0403 
indicates a dominant antiferromagnetic J 2 . Since J 2 cou¬ 
ples triangular units in different c-axis layers, it is re¬ 
sponsible for the TT phase shift between each triangular 
layer. Meanwhile, the nearest neighbor exchange within 
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the triangles may be either ferromagnetic or antiferro¬ 
magnetic. Modaressi, et al analyzed the possible mag¬ 
netic couplings, both direct and indirect, and proposed 
both ferromagnetic and antiferromagnetic contributions 
for Ji and only antiferromagnetic coupling for J2 ^ This 
may allow for a weaker overall Ji compared to J2 (both 
AFM). Furthermore, since the coordination number for 
J2 bonds is 2: = 4 while that for Ji bonds is 2: = 2 , 
comparable AFM Ji and J2 strengths would still pro¬ 
duce a state dominated by J2. For example, they could 
produce a commensurate version of the proposed helical 
state, with all spins aligned on a triangle and a tt phase 
shift from layer to layer (this is symbolically represented 
in Fig. I^b). We used the software package SpinW^ to 
numerically minimize the energy of the AFM Ji — J2 
Heisenberg model for Fe3P0403, starting from random 
spin configurations in a 10 x 10 x 4 supercell with peri¬ 
odic boundary conditions. We found that the aforemen¬ 
tioned “parent” commensurate structure is realized for 
the ratios J2/J1 > 3 . 5 , and a helical incommensurate 
structure like the one refined for Fe3P0403 is realized 
for J2/J1 ^ 3 . The helical modulation therefore arises as 
a compromise for the competing interactions; the com¬ 
mensurate structure internally frustrates the triangular 
units and this can be relieved by slightly canting spins 
away from perfect co-alignment. The pitch of the helical 
modulation, i.e., |(i|, should then be controlled by the ra¬ 
tio of Ji to J2. In Fe3P0403 we found the pitch to be 
86 A; this number can in principle be used to extract 
the ratio of exchange parameters using a more detailed 
calculation. Considering only these completely isotropic 
exchange interactions, the triangular geometry prevents 
this mechanism from producing a preferred modulation 
vector direction in the ab plane. In principle, this could 
lead to domains forming from unique selections from a 
continuous distribution of S directions, which could con¬ 
tribute to the apparent independence of for the mag¬ 
netic refinements on the direction of the modulation vec¬ 
tor, as shown in Fig. However, it is more likely that 
small anisotropies in the interactions, or even entropic 
terms in the free energy, could select preferred directions 
(for instance, leaving a three-fold degeneracy). Such in¬ 
formation would likely be obscured due to powder av¬ 
eraging the neutron diffraction pattern. The simple ar¬ 
guments presented here also do not offer a motivation 
for the modulation vector, to be locked to the helical 
plane direction, n, in contrast to MnSi.^ To first order, 
these two directions (^ and h) may be independent in 
Fe3P0403. 


C. Highly Anisotropic Domains 

Perhaps the most interesting feature of the magnetic 
structure adopted by Fe3P0403 is the needle-like AFM 
domains that maintain small correlation lengths (70 A) 



FIG. 9. Premedial lattice in Fe 3 P 0403 . Unit cells repre¬ 
sented in rhombohedral setting, (a) Top right: schematic of 
iron sites in Fe 3 P 0403 connected by red (Ji) and blue (J 2 ) 
bonds. Bottom right: locations of the centers of triangles 
formed by the Ji bonds, i.e., the simple rhombohedral “pre¬ 
medial” lattice. Left: the relative positions of the premedial 
lattice points (black) with the iron sites (blue), (b) abstrac¬ 
tion of the type of commensurate magnetic structure that the 
helical state in Fe 3 P 0403 is based upon. Each triangular unit 
possesses nearly collinear aligned moments; this direction is 
represented by a single arrow on the premedial lattice, and 
forms a G-type antiferromagnet on a simple rhombohedral 
lattice. 

in the ab plane, but extended correlations (> 900 A) 
along c. The needle-shaped correlation volume implies 
that domain walls cost little energy as long as they are 
“vertical”, i.e., parallel to the hexagonal c axis. 

These anisotropic domains are counterintuitive, given 
that the interactions are three-dimensional. As noted in 
the previous section, the magnetic structure implies that 
J2 is strong and antiferromagnetic. J2 connects the Fe^+ 
triangular units three dimensionally. This can be visu¬ 
alized using the “premedial lattice”, in which a lattice 
point is placed at the center of each triangle.!^ The 
premedial lattice for Fe3P0403 is a simple rhombohe¬ 
dral lattice; i.e., a simple cubic lattice compressed along 
a ( 111 ) direction (Fig. [^a). Each bond between near¬ 
est neighbors on the premedial lattice represents two J2 
bonds in the normal lattice. One can visualize the com¬ 
mensurate version of the magnetic structure realized in 
Fe3P0403 by assigning a single arrow for each trian¬ 
gular subunit (premedial lattice point). The result is a 
G-type antiferromagnet on a distorted simple cubic lat¬ 
tice, a canonical three dimensional antiferromagnet (Fig. 

I^b). 

Given these strong three dimensional couplings, why 
does Fe3P0403 adopt such highly anisotropic correla¬ 
tions? The appearance of vertical domain walls in this 
structure can be partially understood as a consequence of 
the frustrated triangular units that decorate the preme¬ 
dial lattice. The frustration arising from Ji - J2 competi¬ 
tion on the normal lattice allows this simple 3 D structure 
to be fractured by domain walls since they are favored 
by the Ji AFM interaction. Consider a simple domain 
wall consisting of a 180 ° spin reorientation at a plane. 
To investigate the most favorable orientation for such 
a domain wall, we have considered the energy cost for 
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FIG. 10. Understanding domain walls parallel to c in 
Fe 3 P 0403 . (a) Cartoon showing a helical incommensurate 
structure similar to the one found in Fe 3 P 0403 which can be 
generated from competing AFM Ji and J 2 , with J 2 > Ji. 
(b) A simple 180° domain wall passing through the mirror 
planes of the Ji triangular units; spins lying on the domain 
wall are rotated by ~ 90° from the local magnetization direc¬ 
tion. (c) Dehnition of angular parameters for the orientation 
of a domain wall, (d) Simulated energy cost for a domain wall 
with orientation angles 0 and 0 as shown in c), for J 2 /J 1 = 4. 
The lowest energy cost is for “vertical” (parallel to c) domain 
walls passing through local triangular mirror planes = 0, 
0 = 30,90,150°). 


a 180° domain wall in the commensurate parent mag¬ 
netic structure by assuming that when the domain wall 
plane intersects a lattice site, the moment on that site 
lies perpendicular to the (arbitrary) direction chosen for 
the collinear spin structure (Fig.j^a and b). The orien¬ 
tation of the domain wall is parameterized by two angles, 
0 and 0, measured from the hexagonal c and b axes, re¬ 
spectively (Fig. [^c). For J 2 /J 1 = 4, i.e., within the 
parameter range found to produce the parent commen¬ 
surate structure of interest, the energy cost per J2 per 
unit volume is shown in Fig. d) as a function of 0 , 
at different 0 values. In general, the lowest energies are 
achieved at low 0 values, corresponding to “vertical” do¬ 
main walls. Furthermore, the minimum energy is along 
the local mirror planes for the triangular units, as illus¬ 
trated at the top of Fig.j^d). Therefore, for the simple 
domain walls considered here, a threefold degeneracy of 
orientations exists. The intersection of three such domain 
walls could occur at a point, producing a line defect. For 
more realistic domain walls, these special points could 
be topological in nature, similar conceptually to the dis¬ 
crete vortex topological defects found in the hexagonal 
manganites.l^ 

Although helical magnetic structures are common in 
frustrated systems, Fe3P0403 is unusual in that it re¬ 


tains only short range correlations at temperatures 40 
times below the ordering transition. In many cases an 
intermediate incommensurate structure will be entropi- 
cally stabilized above a transition to long range commen¬ 
surate order, a so-called “lock-in transition”. An illustra¬ 
tive example is that of LiNiP04. In this orthorhombic 
system, an intermediate helical incommensurate struc¬ 
ture with short range correlations, similar to that found 
in Fe3P0403, is formed in a small temperature range 
above T/v = 19 K, a t which a commensurate AFM struc¬ 
ture is realized.l^^^ Preceding the transition to commen¬ 
surate order, the incommensurate wavevector (modula¬ 
tion length) shrinks as the temperature is lowered, be¬ 
fore “locking in” to zero at T/v- Higher harmonic peaks 
are observed to grow in as the lock-in transition is ap¬ 
proached, and unusual anisotropic AFM domains are 
observed in this low temperature phase.In contrast, 
the short range incommensurate structure in Fe3P0403 
seems to be stable at all temperatures T < 163 K. There 
is no change in the incommensurate modulation (|^|) 
as a function of temperature, from 160 K down to 4.5 
K, and no higher harmonics have been observed. Fur¬ 
thermore, no increase in the magnetic correlation length 
is observed over this temperature range. Therefore, the 
short range incommensurate structure in Fe3P0403 is 
stabilized, either thermodynamically or perhaps kineti- 
cally (if the domain walls are pinned at lattice defects, 
for example). 

In general, the equilibrium stabilization of AFM do¬ 
main walls is not well-understood. In ferromagnets, do¬ 
main walls are stabilized by demagnetization fields as¬ 
sociated with the net magnetic moment of each domain. 
However, in an ideal antiferromagnet there is no such 
long-range potential to drive fractionalization. A possi¬ 
ble explanation, which so far has been the main route for 
theoretical exploration of this problem for general anti- 
ferro magne ts, is the build up of magnetoelastic strain 
fieldsIn Fe3P0403 such strain fields could be the 
source of the small negative thermal expansion we have 
observed below 100 K (Fig. [^. The extent of magnetoe¬ 
lastic coupling in Fe3P0403 may be important to ex¬ 
plore further in order to explain the stabilization of such 
a high density of AFM domain walls. 

Interest in AFM domains and topological defects is 
rapidly growing due to their importance in exchange 
bias effects and with the new proposals for generation 
of antiferromagnetic Skyrmions. However, they are dif¬ 
ficult to directly investigate experimentally due to the 
net zero moment. Several techniques have been success¬ 
fully employed in the past in order to observe AFM do¬ 
mains or their dynamics, for instance, polarized neutron 
tomography,!^ optical second h armon ic generation,!^ 
and photocorrelation spectroscopy.S^EH Given its unusu¬ 
ally high domain wall density and the potential for gen¬ 
erating topological defects, it would be of great interest 
to explore the domain pattern in Fe3P0403 using these 
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techniques. The development of single crystal samples 
will be of utmost importance for this effort, as well as 
for future neutron scattering studies on this unusual ma¬ 
terial. 


V. CONCLUSIONS 

We have investigated the magnetic and structural 
properties of Fe3P0403 using thermodynamic probes, 
synchrotron X-ray diffraction, and neutron powder 
diffraction. We find that the coexistence of broad, flat- 
topped magnetic peaks with sharp magnetic peaks, as 
observed by neutron powder diffraction, is best described 
by an antiferromagnetic incommensurate helical struc¬ 
ture with highly anisotropic domains that develops below 
T/v = 163 K. The domains are small within the hexag¬ 
onal ab plane, extending to only 70 A, while extending 
to at least 900 A along the hexagonal c-axis, implying 
a high density of ^erticaF domain walls {pDW = 0.65 
nm“^). Stabilization of these antiferromagnetic domain 
walls may occur through magnetoelastic strain; our tem¬ 
perature dependent structure refinements indicate nega¬ 
tive thermal expansion below T/v- 

The magnetic moment determined from our magnetic 
structure refinements, p = 4.14(2)/iB, is consistent with 
the effective moment determined from paramagnetic sus¬ 
ceptibility on the magnetically dilute Fe(3_a,)Gaa,P0403 
series; Peff/Fe ^ b.lpB as the high temperature suscepti¬ 
bility reaches the truly paramagnetic regime. This 16% 
reduced moment compared to the S=b/2 expected for 
Fe^+ could be an indication of covalency (ligand charge 
transfer) in Fe3P0403. 

The type of magnetic structure formed by Fe3P0403 
indicates frustration between the antiferromagnetic ex¬ 
change couplings, Ji and J2, with dominant J2. We 
have shown that despite the three dimensional network 
formed by these strongest interactions in Fe3P0403, the 
presence of frustrated triangular Ji subunits can lead 
to anisotropic needle-like antiferromagnetic domains, as 
observed. The frustrated triangular subunits provide a 
natural location for fracturing the long range magnetic 
structure, and the “path of least resistance” is to prop¬ 
agate these defects along the c-axis. 

Given the helical antiferromagnetic nature of the mag¬ 
netic structure, as well as its propensity to form defects 
(domain walls), Fe3P0403 may be of interest to study in 
the context of topological spin textures such as antiferro¬ 
magnetic Skyrmions.^ Of more fundamental importance 
is the presently poorly-understood mechanism for stabi¬ 
lization of domains in antiferromagnets. In this context, 
Fe3P0403 provides an unusual test case in which a high 
density of small domains form and persist to the lowest 
temperatures, which remains unexplained from either an 
equilibrium or kinetic standpoint. For these reasons, the 
magnetic domain structure in Fe3P0403 certainly de¬ 


serves further investigation, particularly in single crystal 
samples, should they become available. 
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Appendix A: Fe( 3 _a;)Gaa;P 0403 solid solution x-ray 
data 


The solid solution series Fe(3_a.)Gaa^P0403 was inves¬ 
tigated using a Scintag X-ray powder diffractometer (Cu- 
Ka radiation). Select diffraction patterns are shown in 
Fig. along with the refined lattice parameters as a 
function of x. These follow a linear trend, as expected 
for solid solutions. The solubility limit is reached between 
X = 2.S and x = 2.9, where the product phase-separates 
into GaP04 and Ga203 with a small percentage (1.95% 
by mole) of a phase with the Fe(3_a,)Gaa,P0403 struc¬ 
ture type. 
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FIG. 11. X-ray analysis of the Fe( 3 _a,)Gaa;P 0403 series, (a) 
X-ray diffraction patterns (Gu-Ko radiation) of select mem¬ 
bers of the series, showing the progressive increase of scat¬ 
tering angles for equivalent peaks (i.e., decrease in lattice pa¬ 
rameters) with increasing Ga content (a;). Between x = 2.8 
and 2.9, the Fe 3 P 0403 structure type cannot be produced 
as the majority phase, despite many regrindings, (b) Linear 
trends in a vs. x and c vs. x, as expected based on Veg- 
ard’s law (errorbars are smaller than the data points). Blue 
shaded region shows the range over which the desired struc¬ 
ture cannot be produced as the majority phase (solubility 
limit reached). The x = 3.0 values are from a high pressure 
synthesis of Ga 3 P 0403 from Ref. [44l 
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FIG. 12. Gomparison of the peak shape for the lowest angle 
broad feature as observed with the POWGEN instrument 
to two models, (a) The same helical model as proposed for 
Fe 3 P 0403 (Fig.®, but combining a nearly continuous distri¬ 
bution of wavevectors. (b) Square wave model with up to the 
fifth harmonic shown (black, blue and magenta ticks show 
locations of 1^^, 3^^^^ and 5^^ harmonic peaks respectively). 


Appendix B: Further details of neutron powder 
diffraction peak shape 


The peakshape of the lowest Q feature in the POW¬ 
GEN NPD pattern for Fe3P0403 was further investi¬ 
gated using two additional models. The first is the same 
helical model described in the main text, but with sev¬ 
eral choices for the modulation vector superimposed (see 
inset of Fig. 12 for sampling of ((5^, Si), 1.5) wavevectors). 


The result is a nearly continuous distribution of peaks 
within the broad feature (Fig. a). This model still 
requires a Gaussian peak broadening of FWHM = 0.03 
to account for the slope of the edges of the broad 
feature. Note that no magnetic form factor is applied to 
the calculated intensities. 

The second model investigated is an incommensu¬ 
rate square wave model. The square wave would pro¬ 
duce higher harmonics (i.e. kn = {nSa,nSi),1.5), with n 
odd) in the Fourier expansion of the periodic structure, 
with intensities reduced by a factor of 1/n for the 
harmonic peaks. Fig. b) shows how a hypothetical 
square wave structure compares to the measured peak- 
shape uisng ki = (0.015,0.008,1.5), including the 1^^, 
3^^^ and 5^^ harmonics. 


Appendix C: Alternative Conical Magnetic 
Structure for Fe 3 P 0403 

The NPD pattern for Fe3P0403 can be fit by either a 
helical structure or a conical structure. These structures 
are related by the definition of a conical axis direction 
(n) and an opening angle of the cone (/d). The best fits 
have h within the ab plane, and P = 70° (conical) or 
P = 90° (helical). These refinements are shown for h = a 
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in F ig. p!3| and the magnetic structures are illustrated in 
Fig. |14| The helical structure is the same one displayed 
in Fig. and discussed in the text. Note that neither 
structure produces a net moment. Furthermore, in the 
conical structure, the moments connected by J 2 form an 
angle of ^ 140° with one another, while those connected 
by J 2 within the helical structure form an angle of ^ 
175°. Therefore an isotropic AFM J 2 interaction favors 
helical structure over the conical one. 


opening angle (|8) 


Conical {(3 = 70°,h = a) 



cone axis (n) 

Helical {(3 = 90°,h = a) 




Q(A-^) 



FIG. 13. Left: Refinement of a conical structure with open¬ 
ing angle 70° and cone axis along a, which produces a good 
fit to the NPD pattern. Right: Refinement of a helical struc¬ 
ture with opening angle 90° and cone axis along a (the same 
refinement as shown in Fig. [^. 


FIG. 15. Susceptibility vs. temperature in a 10000 Oe field 
for Fe 3 P 0403 from 1.8 K to 900 K. Measurements conducted 
under different settings are shown in different colors. Purple: 
High temperature susceptibility (300 K - 900 K), shifted by 
+0.00011 emu Oe“^molFe“^. Orange: zero-field cooled low 
temperature susceptibility. Blue: field cooled susceptibility. 


a) Conical {(3^70°,h = a) b) Helical (/3=90°, h = £ 



Conical (j0=7O°) 


1 — ^ 

a-#— i 




Helical ()0=9O°) 

77 


H 






FIG. 14. (a) Representations of a conical structure with open¬ 
ing angle 70° and cone axis along a, which produces a good 
fit to the NPD pattern (top: showing only the first layer of 
triangles in the unit cell) (b) Representations of the helical 
structure with opening angle 90° and cone axis along a (the 
same structure as shown in Fig. and discussed in the main 
text.) 


Appendix D: Full y vs. T for Fe 3 P 0403 

The magnetic susceptibility of Fe3P0403, measured 
under an applied field of 1 T, is shown in Figure The 
high temperature (300 K to 900 K) and low temperature 
(1.8 K to 300 K) data are shown in different colors, since 
the measurements were collected in a different sample 
holder. The high temperature data was collected on a 
sample that was fixed to a heater using Zircar cement. 
The low temperature measurement was conducted on 
a sample encapsulated in a plastic holder. To achieve 
quantitative agreement in the overlapping region, the 
high temperature M/H was shifted by +0.00011 emu 
Oe“^mol“^ (5.5%). The lower absolute susceptibility ob¬ 
tained from the high temperature data (as compared to 
the lower temperature data) may be attributable to a 
diamagnetic contribution from the Zircar cement. 


Appendix E: Isothermal magnetization in Fe 3 P 0403 

The magnetization of Fe3P0403 was measured as a 
function of magnetic field at several temperatures from 
300 K down to 2 K. At each temperature the field was 
increased from zero to 9 T, decreased to -9 T, and then 
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A^oH (T) 


FIG. 16. Isothermal magnetization at T = 2 K, 160 K, and 
300 K. A full field-sweep (OT ^ 9T ^ -9T ^ 9T) at each 
temperature is shown. 


increased back to 9 T. All repeated parts of the M vs. 
/ioH curve follow the same line, indicating no observable 
hysteresis at any measured temperature. The M vs. /roH 
curves for three temperatures are shown in Fig. It is 
noteworthy that the curves are linear up to 9 T at all 
temperatures, and only a small fraction of the full mo¬ 
ment is observed (4.2 jib expected at saturation, and ^ 
0.04 /i^ observed at 300 K and 9 T). These observations 
are consistent with very strong antiferromagnetic corre¬ 
lations, as inferred based on the absence of Curie-Weiss 
behavior up to 900 K, as well as the diffuse magnetic 
scattering seen in neutron diffraction at 300 K (see main 
text). 
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